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When a large number TV of independent diffusing particles are placed upon a site of a d- 
dimensional Euclidean lattice randomly occupied by a concentration c of traps, what is the mth 
moment (t™ N ) of the time tj t ?f elapsed until the first j are trapped? An exact answer is given in 
terms of the probability $a/ (t) that no particle of an initial set of M = N, N — 1, . . . ,N — j particles 
is trapped by time t. The Rosenstock approximation is used to evaluate 3>Af(£), and it is found 
that for a large range of trap concentracions the mth moment of tj t N goes as x~ m and its variance 
as x~ 2 , x being ln 2//d (l — c)mJV. A rigorous asymptotic expression (dominant and two corrective 
terms) is given for {t™ N ) for the one-dimensional lattice. 

PACS numbers: 05.40.Fb,66.30.-h 



I. INTRODUCTION 



Statistical problems related to the diffusion of a single random walker in a medium with traps have been subject 
of intense research during the last decades |l|-|8| . Usually it is assumed that the statistical properties of this single 
(N = 1) random walker are representative of the statistical ensemble. However, there are multiparticle (N > 1) 
problems that can not be analyzed in terms of the single walker theory. An example is the number Sn (t) of distinct 
sites visited (or territory explored) up to time t by N independent random walkers all starting from the same origin 
P~p^|. Another multiparticle problem of interest, that as we will see is closely related to that of the territory explored, 
is the description of the order statistic of the diffusion processes, i.e., the estimate of the time at which the jth particle 
of an initial set of N particles all starting from the same origin is trapped. 

The order-statistic problem when the traps are arranged on a (hyper) sphere (i.e., an absorbing boundary at a fixed 
distance) has been thoroughly studied Jl5|-^8[ . In this paper we consider the more difficult problem in which the traps 
are arranged randomly ("the trapping problem") in a d-dimensional Euclidean medium. A related problem, in which 
the A particles are placed on the left of a one-sided random distribution of traps on a one-dimensional lattice, has been 
investigated in Ref. |fT9| . It is interesting to note that recent advances in optical spectroscopy [^0) make it possible 
to monitor this kind of multiparticle dynamic process. Indeed, the simultaneous tracking of N 3> 1 fiuorescently 
labeled particles and the analysis of the diffusive motion of the particles individually is a useful recent technique 
for characterizing heterogeneous microenvironments (in particular for samples dynamically changing in time such as 
biological samples) [|T). An useful feature of the order statistics approach is that it allows one to infer properties 
of the diffusive system (diffusion constant, number of diffusing particles, concentration of traps, effective dimension 
of the diffusive substrate, . . . ) from only the analysis of the behavior of those particles that are first trapped. This 
could be an advantage when it is impractical or impossible to wait until all the reaction is over . 

The order statistics of the trapping process will be described by means of the probability &j,N that j particles 
of the initial set of A diffusing particles have been trapped, and the other N — j have survived, by time t. In 
this paper we consider that all the particles start from the same origin which is free of traps. The moments of the 
time tj^N at which the jth particle of the initial set of N particles is trapped will be calculated from &j,N ■ This 
probability Qj.N will be given in terms of the survival probability $m(0 = ®o.Ai(t) that no particle of an initial set 
of M (M = N, N — 1, . . . ,N — j) has been absorbed by time t. This last quantity will be estimated by means of 
the Rosenstock approximation using expressions for Sjv(i) calculated in Refs. Jl2| Q. It should be noted that our 
approach to the order statistics of the diffusion process in the presence of randomly placed traps is different from that 
used |J - |l8|| for a fixed configuration of traps. What makes the two problems completely different, and hence the way 
of solving them, is that for the case with a given configuration of traps the probability that A particles are trapped 
by time t is simply the Ath power of the probability for a single particle. This simplifying result does not hold when 
many configurations of randomly placed traps are considered. 

The plan of the paper is as follows. In Sec. || we deduce the main formulas that describe the order statistics of 
the trapping process: we relate &j,N to $jv(0 and (t™ N ) to &j^(t). In Sec. Ill we show that the ratio between the 
variance of S^it) and (5jv(*)) 2 goes roughly as (In A) -2 for large A. This suggests that the Rosenstock approximation 
for $Af(0 can lead to good results even when A is large. This is checked in Sec. [V, where we also obtain asymptotic 
expressions (the main term) for (iJV) an0 - the variance. The procedure, based on the Rosenstock approximation, 
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does not provide analytic asymptotic corrective terms for d > 2, although we show that numerical integration is 
feasible leading to excellent results. However, in Sec. [v[ for the one-dimensional lattice we are able to find a 
rigorous asymptotic expression (up to the second-order corrective term) for (i™/v) for large N. Some remarks and the 



conclusions are presented in Sec. VI 



II. ORDER STATISTICS OF THE TRAPPING PROCESS 

Let us first show how to obtain ^j^(t) from $A/(i) with M = N, N — 1, . . . ,N — j. Let ^j.N(t) be the probability 
that j random walkers of the initial set of N have been absorbed by time t by a given configuration of traps and let 
^>(t) = ^o } N(t) be the probability that no single random walker has been absorbed by time t by this configuration of 
traps. Taking into account that (~Y) is the number of different groups of j particles that can be formed from a set of 
N, one finds 



or, using the binomial expansion, 

= f JV ) J2(-iy- m { I )* v '"• (2) 

V J J TO= o \ m / 

Averaging over different configurations of traps and taking into account that 3>jv(£) = ( i & N (t)) and $j t jsr(t) = (\E' Ji Ar(t)), 
we get 

^, N (t) = (-iy(^A^ N (t), (3) 

where the backward difference formula for the jth derivative 

i 



Ai$ N (t) = J2(-ir[')®N- m (t) (4) 

m=0 ^ ' 

has been used. The difference formula in Eq. (^) can be approximated by the derivative: 

when j <C N. Let hj^it) be the probability that the jth absorbed particle of the initial set of N disappears during 
the time interval (t,t + dt]. This quantity is related to Qj^it) by 

d d J 

h j+ i, N (t) = h jtN - j®j,N{t) = —jr Yl $ m,iv(t) (6) 

m=0 

with ho,N = 0. Then, the mth moment of the time at which the jth particle is trapped is given by 

{q N ) = / t m h hN (t)dt (7) 
Jo 

or, using Eq. (^) and integrating by parts, by 

roc 

> = &>+"» / t m -^ jtN (t)dt (8) 

with 

(t? N ) = m / ^-^(tjdfc. (9) 
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Using Eq. (g), Eq. (§) becomes 

(t?+i,N) = + (-^ j (^)^{tT <N )- (10) 

Thus, the order statistics of the trapping problem is described from (t" l N ) only. However, when N and j are large, Eq. 
( |l0| ) is hardly useful numerically because the quantities (t™ N _ r ), that are added and subtracted (and almost cancelled) 
to obtain the jth difference derivative A J (i^uv) nave to be calculated, then, with extraordinary accuracy (which is 
not easy; see Sees. IV and fv]) in order to get a reasonable estimate for the small quantity (— N ) — (t?jy)) 



from the multiplication of the tiny quantity (t™ N ) by the huge binomial coefficient. In Sec. we will show how 
one can surmount, at least partially, this difficulty. 



III. MOMENTS OF THE NUMBER OF DISTINCT SITES VISITED BY N RANDOM WALKERS 

The main purpose of this section is to show that for large N one can approximate (S N (t)) by (SV(t)) 2 . In other 
words, we will show that the ratio Var(Sjv)/ (Sn) 2 is small for large N and that it decreases when N increases. In fact, 
we will show that the simulation results are compatible with the conjecture made in |ll| that [Var(SAr)] 1 ^ 2 / (Sn) ~ 
1/hiN. These results make it very plausible that the Rosenstock approximation of order zero is a reliable method 
for estimating the survival probability $N(t) for not too long times and small concentrations. This will be analyzed 
in Sec. 0. 

The problem of evaluating (S^(t)) for N = 1 has been intensively studied since it was posed by Dvoretzky and 
Erdos p|,p|,p2| . In 1992, Larralde et al. [P,|l0f addressed the problem for N ^> 1 and m = 1 on Euclidean media. They 
disclosed the existence of three time regimes: a very short-time regime [t <C t x ~ ln(iV)/ ln(d)], or regime I, in which 
there are so many particles at every site that all the nearest neighbors of the already visited sites are reached at the 
next step, so that the number of distinct sites visited grows as the volume of an hypersphere of radius t, (Sjv(i)) ~ t d ; 
a very long-time regime (t' x <t), or regime III, that is the final stage in which the walkers move far away from each 
other so that their trails (almost) never overlap and (Sjv(i)) ~ N(Si(t)); and an intermediate regime (f x <C t <C t' x ), 
or regime II, in which there exists a non-negligible probabilitiy of the trails of the particles overlapping. Of course, 
regime III does not exist for d = 1 ( t' x = oo). For d = 2, t' x = exp(TV), and for d > 3, t' x = 7V 2/(d ~ 2) . In the 
simulations carried out in this paper and for the values of N we are interested in (N 3> 1), the particles spend most 
of the time inside regime II, and regime III is never reached. 

For regime II it has been found that |l^Jl3[| 



(S N (t))nS N (t){l-A) (11) 

with 

S N (t) =v (4Dt\nN) d/2 , (12) 

1 oo n 

A = A(N, = In"" N i" 1 In™ In N (13) 



and where, up to second order (n = 2) 
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M • (18) 



Here oj = 7 + In A + fj,ln(d/2), where 7 ~ 0.577215 is the Euler constant, vq is the volume of the hyphersphere with 
unit radius, and A, fi and hi are given in Table | for d = 1, 2 and 3. The diffusion constant is defined by means of 
the Einstein relation 



.,2\ 



2dDt, (19) 
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TABLE I. Parameters that appear in the asymptotic expression of Sjv(t) Eq. (P-ll). The symbol dD 



refers to the d- dimensional simple hypercubic lattice. The parameter p is [2(6D-7r) 3 /3] 1 2 p(0, 1) [^3) where 
p(0,l) = V6/(327r 3 )r(l/24)r(5/24)r(7/24)r(ll/24) ~ 1.516386 §. 
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FIG. 1. {S 2 N )/t 2 versus InN for the two-dimensional lattice when t = 400. The circles are simulation results averaged over 
10 5 configurations for N = 2 2 , . . . , 2 12 and over 10 4 configurations for N = 2 13 , . . . , 2 16 . The lines represent {S N (t)} 2 when the 
main term (dotted line), first-order approximation (dashed line) and second-order approximation (solid line) for (Sjv(t)) are 
used 



for large t, with (r 2 ) being the mean-square displacement of a single random walker. All the numerical results that 
appear in this paper are calculated using D = l/(2d). 

However, the calculation of higher-order moments of S^it) poses a problem of completely different order of mag- 
nitude that still remains unsolved. In Ref. it was conjectured that the functional form of (S^) for Euclidean 



lattices has the same asymptotic structure for all to, namely, the asymptotic structure of Eq. 
conjectured that 



11). Moreover, it was 



Var(5jv) 
(Sn) 2 



1 



In ./V 



1 + 



ln d IniV 
In AT 



(20) 



{Sn} 2 up to 



for large N, where Var(Sjv) = (S%) — (Sn) 2 is the variance of Sjv(i)- Note that Eq. (pOj ) implies (S N ) 
the first-order asymptotic corrective term, as well as Var(Sj\r) ~ t d (ln N) d ~ 2 for large N. 

Simulation data for (S N (t)) for the two-dimensional lattice are compared in Fig. [l] with results obtained from the 
approximation (S% (t)) ~ (SN(t)} 2 in which the zeroth-, first- and second-order asymptotic approximation for (5jv(i)) 
given by Eq. (O) is used. The large difference between the performance of the three asymptotic approximations 
is quite noticeable as well as the excellent result obtained with the second-order approximation. Similar results 
(not shown) are found for d — 1 and d = 3. Figure ^ shows simulation data for the ratio Var(5W) / (Sn) 2 for the 
two-dimensional lattice. We see that for large N this ratio decays roughly as predicted by Eq. (EOT). 



IV. ORDER STATISTICS OF THE TRAPPING PROCESS BY MEANS OF THE ROSENSTOCK 

APPROXIMATION 



The extended Rosenstock approximation (or truncated cumulant expansion) first proposed by Zumofen and Blumen 
is a well-known procedure [0-0 for solving the Rosenstock trapping problem for a single particle (N = 1). Its 
generalization for estimating the (survival) probabilty $ n (t) that no particle of the initial set of N diffusing particles 
has been trapped by time t is straightforward (details can be found in Ref. |^9|) and we will only quote here those 
results that are useful for our objectives. 
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FIG. 2. Simulation results for the ratio {Sn}/ [Var(S'jv)] 1 ^ 2 for the two-dimensional lattice with N = 2 2 ,2 3 ,. 
t = 400. The configurations employed were the same as in Fig. |l[ 



and 



The zeroth-order Rosenstock approximation for estimating $Ar(i) is given by 



*$(t) = e -Ms*(t)> 



(21) 



where A = — ln(l — c) and c is the concentration of traps 

approximation for (Sjsr(t)) [see Eq. (pi)] is used. The first-order Rosenstock approximation is: 



We will write <&$ n \t) to indicate that the nth-order 



(S N (t))lnp[l + 



A Var(SW) 



Then, the error made by using the zeroth-order Rosenstock approximation can be estimated: 

$ N (t) = $<°>(t) [1 + (A 2 Var(SW))] . 



(22) 



(23) 



Thus, the condition A 2 Var(S/v) <C 1 guarantees the good performance of the zeroth-order Rosenstock approximation. 
We have found in Sec. [ill] that Var(5 < Ar) ~ t d (\nN) d ~ 2 so that the zeroth-order Rosenstock approximation works well 
when X 2 t d (\n N) d ~ 2 <C 1. This means that the approximation improves slightly for d = 1 and worsens slightly for 
d = 3 when N increases. For long times, the Rosentock approximation eventually breaks down, the Donsker-Varadhan 
regime settles in, and the survival probability decays in a distinct way known in the literature as Donsker-Varadhan 
behaviour ||. 

Figure p shows the survival probability $j,jv for the two-dimensional lattice obtained from computer simulations 

and from Eq. (||) when the zeroth-order Rosenstock approximation (t) given by Eq. ( |2l"| ) is used. The agreement 
is excellent. 

Now, we evaluate (f^uv) by means of Eq. (|9|) approximating the survival probabity $jv(i) by the the zeroth-order 
Rosenstock approximation (t) for all times : 



(*r: 



JV/ 



exp [-A(Sjv 



(24) 



Notice that with this approximation we are assuming that, in the integration of Eq. (^) that leads to the 
relevant contribution comes from the time interval in which the Rosenstock approximation works. Next, the expression 
for (Sjv(i)) corresponding to the intermediate time regime is used in Eq. (24) for all times. This approximation is 



reasonable if the integrals of mt m 1 $at(^) on the intervals [0, t x ] and [t' x , oo] are negligible versus (t™ N ). Ast x ~ In AT, 



the approximation concerning the first interval is good as long as (lniV) m <C (£iV)j i- e - [ see Eq. ( |25[ ) below], as long 
as A <C (In N)~ d . For t > t' x , one has (Sjsr(t)} ~ N(Si(t)), so that the approximation regarding the interval [t' x , oo] is 
good when Aexp(iV) ^> 1 for d = 2 and AiV 3 ^> 1 for d = 3. Inserting the main asymptotic term of (5iv(t)), namely, 
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t 

FIG. 3. The jih survival probability $,^(4) versus time for (from top to bottom) j = 0,1,2,3,4, with N = 1000 and 
c = 4 x 10~ 4 for the two-dimensional lattice. The lines represent (t), i.e., the zeroth-order Rosenstock approximation with 
(S'jv(i)) given by the second-order asymptotic approximation. The circles are simulation results averaged over 10 6 configurations. 
Inset: <&o,iv(t)- 




InN 

FIG. 4. The function 10 4 /(ti,jv) versus In N for the one-dimensional lattice with concentration of traps c = 8 x 10 3 and 
N = 2 3 , 2 4 , . . . , 2 16 . We plot simulation results averaged over 10 5 configurations(circles) and the asymptotic approximations of 
order (dotted line), order 1 (dashed line) and order 2 (solid line). 



(S N (t)) « v (4:Dt\nN) d / 2 , into Eq. <M) one gets, after a simple integration, a zeroth-order approximation for the 
mth moment of ij. jv : 

JXl + 2m/d) 1 
^VW- ( Awq ) 2 "V<* (4DlniV)m- ^ 

The corrective terms of (S'jv(t)) are not used in Eq. ( pi| ) because their time dependence for the two- and three- 
dimensional cases impedes analytical integration. 

In Figs. (ii,jv) calculated from Eq. ( p5| ) is compared with numerical simulation results. For the two- and three- 
dimensional lattices we also show the results obtained by means of the numerical integration of Eq. (|24|) when the first- 
and second-order asymptotic approximations for (Sjv^)) [cf. Eq. (O) with n — 1 and n — 2, respectively] are used 
for t > t x = (A/ D)h\N . For t <t x , the expression (5jv(t)) = «oi tiorresponding to the short-time regime is used. 
For d = 1, the first- and second-order results are analytical (see Sec. |v|). Figures ^-|| illustrate the great importance 
of the asymptotic corrective terms in the order-statistics quantities. The way in which the lines corresponding to 
the zeroth-order approximation run almost parallel to the simulation results indicates that the corrective term goes 
essentially as (InTV) -1 . This is confirmed in Sec. where it is found that the rigorous asymptotic expression for 
(^jv) f° r the one-dimensional lattice exhibits corrective terms that decay logarithmically with N. 
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FIG. 5. The function 10 3 /(ti,jv) versus In TV for the two-dimensional lattice with c = 4 x 10" 4 and N = 2 3 , 2 4 , . . . , 2 16 . The 
simulation results are averaged over 10 configurations (circles). The dotted line represents the asymptotic approximation of 
order 0. We also plot results obtained by means of the numerical integration of Eq. ( |24| ) when the first-order (dotted line) and 
second-order (solid line) asymptotic approximations for (5jv(t)) are used. 
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FIG. 7. The ratio crj,N/(tj,N) , j = 1 (circles) j = 2 (squares), j ' = 3 (up triangles) j = 4 (down triangles), 
JV = 2 3 ,2 4 ,... ,2 16 , for d = 1 with c = 8 x 10~ 3 (hollow symbols at the top of the figure), d = 2 with c = 4 x 1(T 4 
(filled symbols) and d = 3 with c = 4 x 10~ 5 (symbols with a bar at the bottom of the figure). The simulation results are 
averaged over 10 5 configurations for d — 1 and d = 2 and over 10 4 configurations for d = 3. The lines represent the (main 
order) asymptotic theoretical results, namely, v5 for d = 1, 1 for d = 2 and 0.678968 • ■ • for d = 3. 



From Eq. ( |10| ) and approximating the difference operator A- 7 by the derivative of order j, one finds 

(Au ) 7 (4L>) m J 
for j <C AT, or, in terms of the psi (digamma) function p4[ , 

ft™ ) ~ (t™ ) + m r ( 1 + 2 ™Afl ^)-^(l) f27) 



For 1 <§; j iV one gets 



r(l + 2m/d) 7 + lnj 
'(Au ) 2ro/rf (4D)« (IniV) 1 



because ip(j) — hx(J) + 0(1/ j) and = 7 |24 ]. 

Therefore, the variance ^ N = (t 2 - N ) — (tj,jvv is given by 

2 „ 2 „ r(i + 4/d)-[r(i + 2/d)] 2 

VjN — VlN — Til 9 • \ Z ^> 

3 ' ' (Awo) 4/d (4£>lniV) 2 

Thus, the main asymptotic term of the ratio Cj,N/{tj,N) is independent of j and iV for large N: 

*i.N [r(i + 4/d)-r 2 (i + 2/d)] 1/2 

(t 3 , N ) ~ r(i + 2/d) 1 j 

In Fig. |^ we plot this ratio for the one-, two- and three-dimensional lattices for several values of j and N. The 
simulation results follow closely the theoretical predictions. 

Finally, note that Eqs. (|25|)-(|3(]) are valid for a given d when N — > 00. For a given N and d 00, time regimes 
I and II shrink, i.e., t' x — > 0, so that (Sjv(i)) ~ Aft because (Sift)) ~ < for d > 3. Introducing this relation into Eq. 
(p4|) one gets ftiV) ~ (AN) - " 1 for d — > 00, as expected. 
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V. ORDER STATISTICS OF THE ONE-DIMENSIONAL TRAPPING PROCESS. RIGOROUS RESULTS 



In this section we obtain the order statistics of the trapping process for the one-dimensional lattice from the order 
statistics of the diffusion process in the presence of two fixed traps. Let i™jy( r ) be the TO * n moment of the trapping 
time of the jth particle out of a total of N particles that were initially placed at distance r from a trap in a given 
direction and at a distance greater than r from another trap in the other direction on a line. This quantity is given 

by M 



where t^jv^wi) = Ti^jv(m) + 8j^(m), 



(r) 



4D In kN 



T 3,N{m), 



(31) 



n,Ar(m) = 1 + 



77). 



In kN 
1 



— In In kN — 7 

2 ' 



2 ln^ kN 



1 + 7 + (1 + m) ( — + 7 2 



+ (1 + 777)7! lnln«iV + -(1 + to) In 2 In kN 



O 



ln d In KiV 
In 3 TtiV 



(32) 



8j <N (m) 



777 4— i S n (m) 



In kN ^ 

n—l 



S n (m) = 1 + 



1 



In TdV 



(-l) r 



5n(2) 

(77-1)! 



- In Ih(kJV) 



1 



2(to + 1) 



In 2 hi kN 
In 2 kN 



(33) 



(34) 



and k = 1 / 1/7F . 

In order to get (t™jv), *™jv( r ) ^ s avera g e d over the different positions on which the iV particles can be initially placed 
in an interval free of traps of size L: 



L/2 



dr t™ N (r) 



1 



1 



2m+ 1 KiDlnnN 



2/1! 



— T 



(m). 



(35) 



Next, this quantity is averaged over the size distribution ij(L) = A 2 Lexp[— XL] of the intervals that are free of traps 
P, p. 217] to get the final result 



(t 



3,Nl 



= _ r(l + 2m) T hN {m) 



dL r)(L) q N (L) = 



(36) 



(2A) 2m (4Dln/tJV) m 

for large N and d = 1. In Fig. [|, the theoretical results given by Eq. (36) are compared with simulation data. A 
behaviour very close to that found for traps arranged over a (hyper) spherical surface p6[ is found: the asymptotic 
corrective terms are not at all negligible even for very large values of N, and the second-order asymptotic expression 
is an excellent approximation even for not too large values of N (say, for N > 100). 

Notice that the approximate result obtained in Eq. (||^) agrees, for the one-dimensional case, with the main term of 
Eq. ([56]). This prompts us to investigate to what extent the approximate procedure of Sec. |V| is able to reproduce the 
results of the rigorous asymptotic approach. The answer is that the two approaches lead to the same main term (as 
we have just discovered) and to almost the same first corrective term. For example, using Eq. ([llj) up to first-order 
corrective terms, one gets for j = 1 and m — 1 that 



1 



2A 2 4L> In kN 



In 7r — 27 - 



hi In TV 



41niV 



+ 



(37) 



with a = 0. This expression differs from the rigorous asymptotic formula ( p6| ) in the value of a only: the exact value 



is a = In 2. Finally from Eq. (|36|) one can also obtain for {tj^ 1 N ) — (tJV) the formula ( pq ) which was obtained in 
Sec. |[y] for ci-dimensional media. 

Finally, from Eq. (pq) one gets the variance: 



r(5)W2) 



T 2 (3)r 2 



'j.N 



j.N 



whose main-order asymptotic term reproduces Eq. 



(2A) 4 (4L>lnKA0 2 
91) when d = 1. 



(38) 
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VI. CONCLUSIONS 



The problem addressed in this paper is easy to formulate: When a set of N 3> 1 diffusing particles are placed 
on a site of a rf-dimensional Euclidean lattice ocuppied by a random distribution of static traps, how long is the 
survival time tj t ff of the jth trapped particle? The answer to this order-statistics problem is given in Eq. (^) in terms 
of the probability ^j,N(t) that j particles have been trapped and N — j survive by time t, which, in turn, can be 
expressed [cf. Eq. rt|)] exactly in terms of the survival probabily $m = $o,m that no particle of an initial set of M 
(Al = N, N — 1, . . . ,N — j) has been trapped by time t. 

For the evaluation of 3>jv(*) we resorted to the Rosenstock approximation generalized to the case of N 3> 1 
particles. This approximation is good for small concentrations of traps and small times. Its range of applicabilty 
depends logarithmically on N, improving slightly for d = 1 and worsening slightly for d = 3 when N increases. 
Analytical expressions for the main asymptotic term of mth moment of tj,jsr and its variance crj N for ci-dimensional 
Euclidean media have been found by assuming that the density of traps is such that the contribution of $jv(i) to (t™ N ) 

is negligible in the time regimes I and III. It was found that (i™jv) ~ (A 2 / d hxN) ™ and that the ratio <Tj.N / (tj.N) is 
not at all negligible. In fact <Jj t N is larger than the difference (tj+i^) — (tj,N), which implies that it is not possible 
to infer with certainty the order j of a trapped particle from the time at which it is trapped. However, this ratio 
discriminates clearly the dimension of the Euclidean media in which the particles diffuse. This leads us to consider 
the possibility that this ratio could serve to estimate the dimension of fractal (disordered) media in a dynamical way. 

For the one-dimensional lattice, the previous solution of the order-statistic diffusive problem for a given configuration 
(no randomly distributed) of traps has been used to obtain second-order asymptotic rigourous expressions for {t™ N ) 
and the variance cf^ N . For d > 2 we resorted to numerical integration to obtain higher-order estimates. This 
numerical procedure leads to excellent results, but it is limited to not too large values of N and j because otherwise 
the binomial term that appears in Eq. (||) [or in Eq. (|l^)] becomes intractably large. In all the cases studied, there 
became clear the great importance of the corrective terms in the asymptotic expressions of the moments of the order- 
statistics quantities since the mth corrective term decay mildly as roughly the mth power of the logarithm of N. 
This characteristic behaviour is shared with other cases with different configurations of traps (e.g., fixed traps) and 
substrates (e.g., fractal media). 

We shall finish by mentioning some open problems. First, it would very interesting to estimate the time tw,N by 
which all the particles are eventually absorbed. Notice that the formulae of Sees. IV and |v| are not suitable for this 
purpose as they are valid for estimating tj t N when j <C N only. Also, it would be interesting to describe the order 
statistic of the trapping problem for a trap concentration small enough for the trappingprocess to take place mainly 
inside the Donsker-Varandhan time regime. The recent analysis of Barkema et. al. |j on the crossover from the 
Rosenstock behaviour to the Donsker-Varandhan behaviour should facilitate this task. Finally, it would be desirable 
to extend the results of the present paper to fractal substrates. To this end, the recent results obtained in Ref. |l4| ] 
on the territory explored by a set of random walkers in fractal media should be very useful. 
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